Mitogenomic exploration supports the historical hypothesis of anthropogenic diffusion of a zoonotic parasite Echinococcus multilocularis.

Animal movement across regions owing to human activity can lead to the introduction of pathogens, resulting in disease epidemics with medical and socioeconomic significance. Here, we validated the hypothesis that human activity, such as the transportation of infected animals, has played a significant role in introducing the zoonotic parasite Echinococcus multilocularis into Hokkaido, Japan, by synthesizing and evaluating parasite genetic data in light of historical records. Our analysis indicates that a major genetic group in Hokkaido originated from St. Lawrence Island, USA, which is in accordance with the route suggested by historical descriptions. Moreover, we identified a minor genetic group closely related to parasites found in Sichuan, China. This fact implies that parasite invasion in Japan may result from complex and inadvertent animal translocations. These findings emphasize the anthropogenic impacts on zoonotic parasite spread and provide a crucial perspective for preventing future potential epidemics.

The dynamics of the microbiome in Ixodidae are shaped by tick ontogeny and pathogens in Sarawak, Malaysian Borneo.

Tick-borne diseases have recently been considered a potential emerging public health threat in Malaysia; however, fundamental studies into tick-borne pathogens and microbiome appear limited. In this study, six tick species (Ixodes granulatus, Haemaphysalis hystricis, Haemaphysalis shimoga, Dermacentor compactus, Dermacentor steini and Dermacentor atrosignatus) collected from two primary forests and an oil palm plantation in Sarawak, Malaysian Borneo, were used for microbiome analysis targeting bacterial 16S rDNA using next-generation sequencing (NGS). In addition, bacterial species were further characterized in conventional PCRs to identify potential pathogens. Sequences generated from NGS were first filtered with the Decontam package in R before subsequent microbial diversity analyses. Alpha and beta analyses revealed that the genus Dermacentor had the highest microbial diversity, and H. shimoga significantly differed in microbial composition from other tick species. Alpha and beta diversities were also significantly different between developmental stages of H. shimoga. Furthermore, we observed that some bacterial groups were significantly more abundant in certain tick species and developmental stages of H. shimoga. We tested the relative abundances using pairwise linear discriminant analysis effect size (LEfSe), which also revealed significant microbial composition differences between Borrelia-positive and Borrelia-negative I. granulatus ticks. Finally, pathogenic and potentially pathogenic bacteria circulating in different tick species, such as Rickettsia heilongjiangensis, Ehrlichia sp., Anaplasma sp. and Bartonella spp. were characterized by PCR and sequencing. Moreover, Coxiella and Francisella-like potential symbionts were identified from H. shimoga and D. steini, respectively. More studies are required to unravel the factors associated with the variations observed in this study.

Dietary niche breadth influences the effects of urbanization on the gut microbiota of sympatric rodents.

Cities are among the most extreme forms of anthropogenic ecosystem modification, and urbanization processes exert profound effects on animal populations through multiple ecological pathways. Increased access to human-associated food items may alter species' foraging behavior and diet, in turn modifying the normal microbial community of the gastrointestinal tract (GIT), ultimately impacting their health. It is crucial we understand the role of dietary niche breadth and the resulting shift in the gut microbiota as urban animals navigate novel dietary resources. We combined stable isotope analysis of hair and microbiome analysis of four gut regions across the GIT to investigate the effects of urbanization on the diet and gut microbiota of two sympatric species of rodents with different dietary niches: the omnivorous large Japanese field mouse (Apodemus speciosus) and the relatively more herbivorous gray red-backed vole (Myodes rufocanus). Both species exhibited an expanded dietary niche width within the urban areas potentially attributable to novel anthropogenic foods and altered resource availability. We detected a dietary shift in which urban A. speciosus consumed more terrestrial animal protein and M. rufocanus more plant leaves and stems. Such changes in resource use may be associated with an altered gut microbial community structure. There was an increased abundance of the presumably probiotic Lactobacillus in the small intestine of urban A. speciosus and potentially pathogenic Helicobacter in the colon of M. rufocanus. Together, these results suggest that even taxonomically similar species may exhibit divergent responses to urbanization with consequences for the gut microbiota and broader ecological interactions.

Comparative mitogenomics elucidates the population genetic structure of Amblyomma testudinarium in Japan and a closely related Amblyomma species in Myanmar.

Ticks are the second most important vector capable of transmitting diseases affecting the health of both humans and animals. Amblyomma testudinarium Koch 1844 (Acari: Ixodidae), is a hard tick species having a wide geographic distribution in Asia. In this study, we analyzed the composition of A. testudinarium whole mitogenomes from various geographical regions in Japan and investigated the population structure, demographic patterns, and phylogeographic relationship with other ixodid species. In addition, we characterized a potentially novel tick species closely related to A. testudinarium from Myanmar. Phylogeographic inference and evolutionary dynamics based on the 15 mitochondrial coding genes supported that A. testudinarium population in Japan is resolved into a star-like haplogroup and suggested a distinct population structure of A. testudinarium from Amami island in Kyushu region. Correlation analysis using Mantel test statistics showed that no significant correlation was observed between the genetic and geographic distances calculated between the A. testudinarium population from different localities in Japan. Finally, demographic analyses, including mismatch analysis and Tajima's D test, suggested a possibility of recent population expansion occurred within Japanese haplogroup after a bottleneck event. Although A. testudinarium has been considered widespread and common in East and Southeast Asia, the current study suggested that potentially several cryptic Amblyomma spp. closely related to A. testudinarium are present in Asia.

Novel symbionts and potential human pathogens excavated from argasid tick microbiomes that are shaped by dual or single symbiosis.

Research on vector-associated microbiomes has been expanding due to increasing emergence of vector-borne pathogens and awareness of the importance of symbionts in the vector physiology. However, little is known about microbiomes of argasid (or soft-bodied) ticks due to limited access to specimens. We collected four argasid species (Argas japonicus, Carios vespertilionis, Ornithodoros capensis, and Ornithodoros sawaii) from the nests or burrows of their vertebrate hosts. One laboratory-reared argasid species (Ornithodoros moubata) was also included. Attempts were then made to isolate and characterize potential symbionts/pathogens using arthropod cell lines. Microbial community structure was distinct for each tick species. Coxiella was detected as the predominant symbiont in four tick species where dual symbiosis between Coxiella and Rickettsia or Coxiella and Francisella was observed in C. vespertilionis and O. moubata, respectively. Of note, A. japonicus lacked Coxiella and instead had Occidentia massiliensis and Thiotrichales as alternative symbionts. Our study found strong correlation between tick species and life stage. We successfully isolated Oc. massiliensis and characterized potential pathogens of genera Ehrlichia and Borrelia. The results suggest that there is no consistent trend of microbiomes in relation to tick life stage that fit all tick species and that the final interpretation should be related to the balance between environmental bacterial exposure and endosymbiont ecology. Nevertheless, our findings provide insights on the ecology of tick microbiomes and basis for future investigations on the capacity of argasid ticks to carry novel pathogens with public health importance.

High infection rate of tick-borne protozoan and rickettsial pathogens of cattle in Malawi and the development of a multiplex PCR for Babesia and Theileria species identification.

Malawi has an estimated cattle population of 1,884,803 heads, the indigenous Malawi zebu breed accounts for 91.2%, while the exotic and crossbred accounts for the remaining 8.8%. Although ticks and tick-borne diseases are widespread in Malawi, no molecular study has been conducted to investigate the tick-borne Anaplasmataceae and piroplasms infecting cattle. To provide an insight into the current status of tick-borne pathogens (TBPs) of cattle, a molecular survey was conducted in the central and southern regions of Malawi. A total of 191 cattle of which 132 were Malawi zebu, 44 were Holstein Friesian and 15 were Holstein-Friesian/ Malawi zebu crosses were screened for Anaplasmataceae and piroplasms using the heat shock protein groEL gene and 18S rDNA, respectively. A new 18S rDNA multiplex PCR assay was designed for Babesia and Theileria species identification without sequencing. Overall, 92.3% (n = 177) of the examined animals were infected with at least one TBP. Anaplasmataceae-positive rate was 57.6% (n = 110) while for piroplasms it was 80.1% (n = 153). The detected Anaplasmataceae were Anaplasma bovis 2.6% (n = 5), Anaplasma marginale 24.6% (n = 47), Anaplasma platys-like 13.6% (n = 26), uncharacterized Anaplasma sp. 14.1% (n = 27), and uncharacterized Ehrlichia sp. 16.2% (n = 31). The detected piroplasms were Babesia bigemina 2.6% (n = 5), Theileria mutans 73.8% (n = 141), Theileria parva 33.0% (n = 63), Theileria taurotragi 12.6% (n = 24), and Theileria velifera 53.4% (n = 102). Mixed infection rate was found in 79.6% (n = 152) of the samples analyzed. This study has shown a high burden of TBPs among cattle in Malawi which highlights the need to conceive new methods to control ticks and TBPs in order to improve animal health and productivity. The newly developed multiplex PCR assay would be a useful tool especially in resource limited settings where sequencing is not available and when mixed infections are expected.

Reconstruction of mitochondrial genomes from raw sequencing data provides insights on the phylogeny of Ixodes ticks and cautions for species misidentification.

High-throughput sequencing (HTS) technology has profoundly been involved in sequencing whole genomes of several organisms in a fast and cost-effective manner. Although HTS provides an alternative biomonitoring method to the time-consuming and taxonomy-expertise dependent morphological approach, still we cannot rule out the possibility of the impediment and misidentification biases. In this article we aim to retrieve whole mitochondrial genome (mitogenome) sequences from publicly available raw sequencing data for phylogenetic comparison of Ixodes persulcatus. For this comparison, we sequenced whole mitogenomes of four I. persulcatus ticks from Japan and constructed mitogenomes from raw sequencing data of 74 I. persulcatus ticks from China. Bayesian phylogenetic trees were inferred by the concatenated fifteen mitochondrial genes. We further tested our results by the phylogenetic analysis of cytochrome c oxidase subunit 1 (cox1) gene and internal transcribed spacer 2 (ITS2) sequences. Our findings showed that 70 constructed mitogenomes from China were clustered with the sequenced four mitogenomes of I. persulcatus from Japan. We also revealed that mitogenome sequences retrieved from two data sets CRR142297 and CRR142298 were clustered with Ixodes nipponensis. Moreover, other two mitogenome sequences from CRR142310 and CRR142311 formed a clade with Ixodes pavlovskyi. The phylogenetic analysis of cox1 gene and ITS2 sequences confirmed the identification errors of these four samples. The overall phylogenetics in our study concluded that accurate morphological identification is necessary before implementing HTS to avoid any misidentification biases.

Comparing the gut microbiome along the gastrointestinal tract of three sympatric species of wild rodents.

Host-microbe interactions within the gastrointestinal tract (GIT) play a pivotal role in shaping host physiology, ecology, and life history. However, these interactions vary across gut regions due to changes in the physical environment or host immune system activity, thereby altering the microbial community. Each animal species may harbor their own unique microbial community due to host species-specific ecological traits such as dietary habits, micro-habitat preferences, and mating behavior as well as physiological traits. While the gut microbiota in wild animals has received much attention over the last decade, most studies comparing closely related species only utilized fecal or colon samples. In this study, we first compared the gut microbial community from the small intestine, cecum, colon, and rectum within three sympatric species of wild rodents (i.e. Apodemus speciosus, A. argenteus, and Myodes rufocanus). We then compared each gut region among host species to determine the effect of both gut region and host species on the gut microbiota. We found that the small intestine harbored a unique microbiome as compared to the lower GIT in all three host species, with the genus Lactobacillus in particular having higher abundance in the small intestine of all three host species. There were clear interspecific differences in the microbiome within all gut regions, although some similarity in alpha diversity and community structure within the small intestine was found. Finally, fecal samples may be appropriate for studying the lower GIT in these species, but not the small intestine.

Applications of Blocker Nucleic Acids and Non-Metazoan PCR Improves the Discovery of the Eukaryotic Microbiome in Ticks.

Ticks serve as important vectors of a variety of pathogens. Recently, the viral and prokaryotic microbiomes in ticks have been explored using next-generation sequencing to understand the physiology of ticks and their interactions with pathogens. However, analyses of eukaryotic communities in ticks are limited, owing to the lack of suitable methods. In this study, we developed new methods to selectively amplify microeukaryote genes in tick-derived DNA by blocking the amplification of the 18S rRNA gene of ticks using artificial nucleic acids: peptide nucleic acids (PNAs) and locked nucleic acids (LNAs). In addition, another PCR using non-metazoan primers, referred to as UNonMet-PCR, was performed for comparison. We performed each PCR using tick-derived DNA and sequenced the amplicons using the Illumina MiSeq platform. Almost all sequences obtained by conventional PCR were derived from ticks, whereas the proportion of microeukaryotic reads and alpha diversity increased upon using the newly developed method. Additionally, the PNA- or LNA-based methods were suitable for paneukaryotic analyses, whereas the UNonMet-PCR method was particularly sensitive to fungi. The newly described methods enable analyses of the eukaryotic microbiome in ticks. We expect the application of these methods to improve our understanding of the tick microbiome.

Exploring Prokaryotic and Eukaryotic Microbiomes Helps in Detecting Tick-Borne Infectious Agents in the Blood of Camels.

Dromedary camels (Camelus dromedarius) are widely distributed in Africa, the Middle East and northern India. In this study, we aimed to detect tick-borne pathogens through investigating prokaryotic and eukaryotic microorganisms in camel blood based on a metagenomic approach and then to characterize potentially pathogenic organisms using traditional molecular techniques. We showed that the bacteria circulating in the blood of camels is dominated by Proteobacteria, Bacteroidetes, Firmicutes and Actinobacteria. At the genus level, Sediminibacterium, Hydrotalea, Bradyrhizobium and Anaplasma were the most abundant taxa. Eukaryotic profile was dominated by Fungi, Charophyta and Apicomplexa. At the genus level, Theileria was detected in 10 out of 18 samples, while Sarcocystis, Hoplorhynchus and Stylocephalus were detected in one sample each. Our metagenomic approach was successful in the detection of several pathogens or potential pathogens including Anaplasma sp., Theileria ovis, Th. separata, Th. annulate, Th. mutans-like and uncharacterized Theileria sp. For further characterization, we provided the partial sequences of citrate synthase (gltA) and heat-shock protein (groEL) genes of Candidatus Anaplasma camelii. We also detected Trypanosoma evansi type A using polymerase chain reaction (PCR) targeting the internal transcribed spacer 1 (ITS1) region. This combined metagenomic and traditional approach will contribute to a better understanding of the epidemiology of pathogens including tick-borne bacteria and protozoa in animals.

Spiroplasma Infection among Ixodid Ticks Exhibits Species Dependence and Suggests a Vertical Pattern of Transmission.

Members of the genus Spiroplasma are Gram-positive bacteria without cell walls. Some Spiroplasma species can cause disease in arthropods such as bees, whereas others provide their host with resistance to pathogens. Ticks also harbour Spiroplasma, but their role has not been elucidated yet. Here, the infection status and genetic diversity of Spiroplasma in ticks were investigated using samples collected from different geographic regions in Japan. A total of 712 ticks were tested for Spiroplasma infection by PCR targeting 16S rDNA, and Spiroplasma species were genetically characterized based on 16S rDNA, ITS, dnaA, and rpoB gene sequences. A total of 109 samples originating from eight tick species were positive for Spiroplasma infection, with infection rates ranging from 0% to 84% depending on the species. A linear mixed model indicated that tick species was the primary factor associated with Spiroplasma infection. Moreover, certain Spiroplasma alleles that are highly adapted to specific tick species may explain the high infection rates in Ixodes ovatus and Haemaphysalis kitaokai. A comparison of the alleles obtained suggests that horizontal transmission between tick species may not be a frequent event. These findings provide clues to understand the transmission cycle of Spiroplasma species in wild tick populations and their roles in host ticks.

Anthropogenic interferences lead to gut microbiome dysbiosis in Asian elephants and may alter adaptation processes to surrounding environments.

Human activities interfere with wild animals and lead to the loss of many animal populations. Therefore, efforts have been made to understand how wildlife can rebound from anthropogenic disturbances. An essential mechanism to adapt to environmental and social changes is the fluctuations in the host gut microbiome. Here we give a comprehensive description of anthropogenically induced microbiome alterations in Asian elephants (n = 30). We detected gut microbial changes due to overseas translocation, captivity and deworming. We found that microbes belonging to Planococcaceae had the highest contribution in the microbiome alterations after translocation, while Clostridiaceae, Spirochaetaceae and Bacteroidia were the most affected after captivity. However, deworming significantly changed the abundance of Flavobacteriaceae, Sphingobacteriaceae, Xanthomonadaceae, Weeksellaceae and Burkholderiaceae. These findings may provide fundamental ideas to help guide the preservation tactics and probiotic replacement therapies of a dysbiosed gut microbiome in Asian elephants. More generally, these results show the severity of anthropogenic activities at the level of gut microbiome, altering the adaptation processes to new environments and the subsequent capability to maintain normal physiological processes in animals.

Amblyomma testudinarium infestation on a brown bear (Ursus arctos yesoensis) captured in Hokkaido, a northern island of Japan.

The tick Amblyomma testudinarium Koch, 1844 (Acari: Ixodidae) is known as a vector of several pathogens such as Rickettsia tamurae and severe fever with thrombocytopenia syndrome (SFTS) virus. This tick species is present in many Asian countries, including Japan, where its distribution is limited to the warm areas of Kanto region and the southwestern region. The present study reports the recovery of a partially engorged A. testudinarium from a wild brown bear captured in Shari town, Hokkaido. In addition to morphological identification, the specimen was genetically characterized by the complete mitochondrial genome sequencing. The results showed that the length of the obtained mitogenome is 14,835 bp that encodes 13 protein-coding, two ribosomal RNA (rRNA) (12S and 16S), and 22 transfer RNA genes with two non-coding control regions. The phylogenetic analysis indicated that our sample clustered with A. testudinarium from Nara, Japan, but separated from A. testudinarium from China. Although the introduction of the tick through livestock transportation cannot be ruled out, the detection of A. testudinarium in Hokkaido prefecture, which is separated from the main island where A. testudinarium is present in the south, may suggest the introduction by migratory birds. This study provides important insights on the distribution and host range of A. testudinarium. This will be useful for the future taxonomic analysis of ticks based on the complete mitogenome sequencing. To our knowledge, this is the northernmost detection point of the tropical tick A. testudinarium.

Molecular identification and characterization of piroplasm species in Hokkaido sika deer (Cervus nippon yesoensis), Japan.

Babesia and Theileria species are tick-borne protozoan parasites that have a veterinary and zoonotic importance. In order to investigate the prevalence and genetic diversity of these parasites, a total of 269 sika deer blood DNA samples collected from Hokkaido, Japan, were examined for Babesia and Theileria species by touch-down PCR targeting the 18S rRNA gene. Reverse line blot (RLB) hybridization was then used to detect 12 piroplasm species. The results revealed that 95.5% (257/269), 94.1% (253/269), 14.1% (38/269), 87.7% (236/269) and 11.5% (31/269) of the examined PCR products hybridized with the probes which were designed to detect all Babesia and Theileria spp., all Theileria spp., all Babesia spp., Theileria sp. Thrivae and Babesia divergens-like, respectively. The 18S rRNA gene partial sequences were divided into Theileria sp. Thrivae, T. capreoli, B. divergens-like and an undescribed Babesia species. This study showed the first detection of the undescribed Babesia sp. from Japan. Therefore, more studies are required to understand the ecology of the newly detected tick-borne pathogens in Hokkaido.